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Abstract. The effects of background line opacity (line-blocking) in statistical equilibrium calculations for Fe in late-type stellar 
atmospheres have been investigated using an extensive and up-to-date model atom with radiative data primarily from the iron 
Project. The background metal line opacities have been computed using data from the marcs stellar model atmospheres. While 
accounting for this line opacity is important at solar metallicity, the differences between calculations including and excluding 
line-blocking at low metallicity are insignificant for the non-local thermodynamic equilibrium (non-LTE) abundance corrections 
for Fe i lines. The line-blocking has no impact on the non-LTE effects of Fe n lines. The dominant uncertainty in Fe non-LTE 
calculations for metal-poor stars is still the treatment of the inelastic H i collisions, which have here been included using scaling 
factors to the classical Drawin formalism, and whether or not thermalisation of the high Fe i levels to Fe n ground state should 
be enforced. Without such thermalisation, the Fei non-LTE abundance corrections are substantial in metal-poor stars: about 
0.3 dex with efficient (i.e. Drawin-like) Hi collisions and S 0.5 dex without. Without both thermalisation and Hi collisions, 
even Fen lines show significant non-LTE effects in such stars. 
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1. Introduction 

Accurate knowledge of stellar Fe abundances is crucial for our 
understanding of stellar evolution as well as for the study of 
cosmic chemical evolution. Indeed the abundance of iron is 
normally taken as a proxy for the stellar metal content. Recent 
work has been devoted to quantify possible departures from 
local thermodynamic equilibrium (LTE) in the formation of 
Fe I spectral lines but without reaching any consensus to the 
expected magnitude of the non-LTE effects in late-type stars 
(e.g. Gratton et al. lT^l Thevenin & Idiart [T333l Shchukina & 
Trujillo Bueno l200Tl Gratton et al l200Tl Gehren et al. l2001al 
2001b i Korn et al. 2003 Shchuki na et al . 120051 see discussion 
and further references in Asplund 2005). The different results 
are particularly striking at low metallicity where some authors 
find very large non-LTE abundance corrections (> +0.3 dex) 
while others find negligible effects (< +0.1 dex). It is embar- 
rassing that such large discrepancies still exist for the most 
commonly studied element. It certainly raises concern over the 
astronomical interpretations not only for Fe but for other ele- 
ments as well. 
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Since the main non-LTE physical mechanism for Fe - over- 
ionisation - feeds on the UV radiation field it is paramount to 
include the effects of line-blocking in the computations. This 
study has been motivated by the fact that many of the above- 
mentioned non-LTE calculations have been performed with- 
out properly including this line-blocking. For example, back- 
ground line opacities were not accounted for in detail in the 
influential non-LTE study of Thevenin & Idiart (1999). Instead 
artificial multiplication factors were applied to the continuous 
opacities in their calculations in order to simulate the line haze 
in the UV. Could the large non-LTE abundance corrections they 
found be significantly over-estimated because of this shortcom- 
ing? 

Here we attempt to quantify the effect of line-blocking on 
the determination of non-LTE abundance corrections for Fe i 
lines. In order to achieve this, we have modified the widely 
used statistical equilibrium code multi (Carlsson 1986) and as 
such our work will also be of interest for future non-LTE line 
formation studies for other elements. In view of the remaining 
uncertainties attached to the treatment of inelastic H i collisions 
(e.g. Asplund|2005 1, our focus in this initial investigation is not 
on the exact values of the non-LTE abundance corrections but 
on the differential impact of line-blocking in late-type stars of 
different metallicity. We defer a comprehensive investigation 
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of the non-LTE line formation for a large grid of stellar model 
atmospheres to a future study. We also intend to return to the 
problem of Fe line formation using 3D model atmospheres. 



2. Methods 

2.1. The model atom 

2.1 .1 . Energy levels and radiative data 

For the non-LTE Fe i spectral line formation calculations in this 
paper we adopted an updated version of the Fe model atom by 
Thevenin & Idiart ( 1999 1. At the present time the Fe i and Fe n 
systems in the National Institute of Standards and Technology 
(NIST) database 1 consist of 493 and 578 identified levels re- 
spectively (Sugar & Corliss 1985 1. The laboratory analysis by 
Nave et al. ( 1994 1 has further extended the information on the 
Fe i system by bringing the number of identified energy lev- 
els to 846 and the number of identified lines to nearly 9 800. 
In recent years moreover the ongoing international collabora- 
tion known as the ferrum Project (Li et al. 1993 Johansson et 
al. 2002 1 has aimed at expanding and improving the database 
of oscillator strengths and transition probabilities of Fe n lines 
of astrophysical significance by means of both laboratory mea- 
surements and theoretical calculations. However, as the ferrum 
Project has only addressed a selection of Fen bound-bound 
transitions at this stage and is not yet complete we relied on 
the data provided by NIST to construct the Fe n sub-system in 
our work. In practice data storage and computing time impose 
limitations to the complexity of the actual model atom chosen 
for the calculations. We therefore restricted our model to 334 
levels of Fe i from the NIST database up to an excitation poten- 
tial of 6.91 eV (about 1 eV below the first ionisation edge), 189 
levels of Fe n up to 16.5 eV above the Fe i ground state, and the 
ground level of Fe m. In total 3 466 bound-bound radiative tran- 
sitions between Fe i levels and 3 440 between Fe n levels were 
considered, together with the photo-ionisation cross-sections 
for 523 bound-free transitions. 

Spectral line profiles were sampled for typically 70 wave- 
length points and approximated with the standard Voigt func- 
tion. The ^/-values for Fe i lines were taken from Nave et al. 



(19941 and Kurucz & Bell (I1995> databases, while for the os- 
cillator strengths of Fe n lines we referred to the same sources 
as Thevenin & Idiart ( 1999 1. The radiative data for the photo- 
ionisation cross-sections came from the iron Project (Bautista 
I1997> . The detailed bound-free cross-sections were smoothed 
whenever possible to reduce the number of wavelength points 
in the calculations and therefore decrease the computing time. 
When strong resonances were present near the photo-ionisation 
edges we kept the regular frequency sampling. We tabulated 
the cross-sections for bound-free transitions from Fe i levels 
for typically about 100 wavelength points on average up to 300 
wavelength points when high-resolution was necessary. Cross- 
sections for photo-ionisations from Fe n levels were tabulated 
instead for about 30 wavelength points on average. 



Table 1. Stellar parameters for which the marcs model at- 
mospheres used in this paper were computed. The iron abun- 
dance is expressed in the notation [Fe/H] = log(n Fe /n H ) _ 
log(n Fe /n H )©- A solar iron abundance loge Fe Q = 7.50 is as- 
sumed. 



Object 


7W[K] 




[Fe/HJuE 


Sun 


5 780 


4.44 


0.0 


Procyon 


6 500 


4.00 


0.0 


HD 140283 


5 700 


3.70 


-2.5 


G64-12 


6400 


4.10 


-3.3 



2.1.2. Collisional data 



0000 



? 8000 



CD 

CL 



6000 



4000 



SUN VAL-3C 
SUN MARCS 
SUN HM 



-6 



-4 -2 
Log(T 5000 ) 



CD 

1 

O 

CD 
CL 



:oooo 



8000 



6000 



4000 



G64-12 
HD 140283 
Procyon 



-6 



-4 -2 
Log(T 5000 ) 



Fig. 1. Temperature stratification versus standard optical depth 
at A = 5000 A for the three models solar atmosphere adopted in 
this paper (upper panel) and for the marcs model atmospheres 
of Procyon, HD 140283 and G64-12 (lower panel). 



All levels in our Fe model atom were coupled with one 
another via collisional excitation and ionisation by electrons 
and by neutral hydrogen atoms, as these types of collisions 
are believed to be the most relevant in the atmospheres of 
cool stars. Collisional excitations by electrons were imple- 
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mented using van Regemorter's formula ( 1962 > while cross 
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sections for ionisation by electron impact were computed fol- 
lowing Cox (2000 1. All oscillator strengths / were set to a 
minimum value of 10~ 3 in the formulas for bound-bound col- 
lisions with electrons and hydrogen atoms. For the latter type 
of collisions we adopted the approximate description proposed 
by Drawin (1968 1969 1, as done by Steenbock & Holweger 
( 1984 1. However, both laboratory measurements and quantum- 
mechanical calculations for simple atoms indicate that the 
Drawin's recipe fails to reproduce the right order of magni- 
tude for the cross-sections of Hi collisions: in the cases for 
which comparison with experimental data or theoretical cal- 
culations is possible, cross-sections computed with the aid of 
Drawin's formula appear to be overestimated by factors of 
10 - 10 6 depending on the species and transition (Fleck et al. 
fT9l)9l Lambert IT9l%l Belyaev et al. [T9^9l Belyaev & Barklem 
[T9991 Barklem et al. l2003l . 

On the other hand, some recent stellar spectroscopic stud- 
ies suggest that the cross-sections of Hi collisions with Fe 
atoms might actually be reasonably well described by Drawin's 
recipe. Korn et al. (2003 1 found that ionisation equilibrium be- 
tween Fe i and Fe n in cool metal-poor stars is matched consis- 
tently with surface gravities derived from hipparcos parallaxes 
when the Drawinian collisional cross-sections are scaled by a 
factor S h ~ 3 using their particular model dependent r e (j scale 
(based on H lines). Gratton et al. (1999 1 calibrated Hi colli- 
sions with Fe by using RR Lyrae variables at their minimum 
light and adopted an even larger scaling factor Sh ~ 30. The 
observational evidence is, however, not unambiguous as some 
studies have found the presence of large differences between 
Fei and Fen based LTE abundances in halo stars (e.g. Nissen 
et al. ll997l Fuhrmann 179981 Allende Prieto et al. [T999l which 
might indicate that the Hi collisions be inefficient compared 
with Drawin's recipe. Other observational studies suggest van- 
ishing or very small differences between the Fe i and Fe n based 
abundances (e.g. Gratton et al. l^OTTTl Kraft & Ivans l^UOll . The 
efficiency of H i collisions is clearly an open issue which re- 
quires further work both observationally and, most importantly, 
through detailed atomic physics calculations. 

As a reliable and comprehensive description of neutral hy- 
drogen collisions with iron atoms is not available at this stage, 
we treated the efficiency of H i collisions a la Drawin as a free 
parameter throughout our study. In particular, considering the 
large uncertainty on the actual cross-sections of H i collisions, 
we decided to apply multiplication factors Sh = 0.001 and 
Sh = 1 to Drawin's formula and compared the resulting effects 
on the non-LTE calculations in the two cases. This approach is 
justified since the main focus of the present paper is to investi- 
gate the effects of line-blocking rather than providing the final 
word on the size of the Fe non-LTE abundance corrections for 
individual stars. 

2.2. The model atmospheres 

As the purpose of our study was to investigate the effects of 
line-blocking on the non-LTE formation of Fe i spectral lines, 
we have not performed non-LTE calculations for a large grid 
of models. Instead we have selected a sample of four cool stars 



(the Sun, Procyon, HD 140283 and G64-12), covering the most 
relevant parameter space and performed the necessary test cal- 
culations. 

The analysis was carried out using marcs line-blanketed, 
plane-parallel model atmospheres (Asplund et al. 1997) of the 
four objects, computed for the set of stellar parameters indi- 
cated in Table^ an d a micro-turbulence £ = 1.0 kms -1 . We 
kept these effective temperatures and surface gravities constant 
throughout the analysis, as the main purpose of our investiga- 
tion is to study differentially the effect of the inclusion and ex- 
clusion of line-blocking on the derivation of non-LTE correc- 
tions to the iron abundance with respect to the LTE case and 
not to fine-tune stellar parameters. We note that our conclu- 
sions are immune to reasonable changes of the adopted stellar 
parameters. In particular the computed Fe non-LTE abundance 
corrections are affected by about 0.01 dex or less when a micro- 
turbulence £ = 2.0 kms -1 is used instead. 

Finally, we adopted for the Sun the solar chemical com- 
position by Grevesse & Sauval (1998 1, i.e. assuming an iron 
abundance of log e Fe = 7.50 2 . The chemical composition of 
the other models was scaled with respect to solar according to 
their LTE iron abundance and assuming an or-enhancement of 
[a/Fe] = 0.4 dex for the two metal-poor stars in the sample. 

Along with the marcs solar atmosphere, we considered two 
additional models for the Sun: the semi-empirical Holweger- 
Miiller (hm) model (Holweger & Miiller H"974> and the quiet- 
Sun val-3c model by Vernazza et al. dl981> . In this way we in- 
tended to cover - at least for the Sun - typical model-to-model 
variations in the analysis. The temperature stratifications for 
the three solar models are shown in Fig.^together with mod- 
els for the other three stars. The marcs and hm models have 
similar stratifications although the latter has slightly warmer 
temperatures in the line forming regions. The non-LTE val-3c 
model departs from the other two solar models mainly by pre- 
senting a chromospheric temperature rise in its upper layers 
and reaching significantly lower temperatures at optical depths 
-4 S5 log(T5ooo) ^ -1.5. For a quantitative description of how 
structural differences arise between LTE and non-LTE mod- 
els of the solar photosphere we refer to the work by Rutten & 
Kostik JT9821 and Rutten (IT9881 

2.3. Radiative transfer calculations 

We computed the emerging Fe spectrum for each of the model 
atmospheres by solving the statistical equilibrium and radia- 
tive transfer equations with the non-LTE code multi (Carlsson 
1986 1, version 2.2. We introduced however some modifications 
concerning the treatment of background opacities, multi's stan- 
dard background opacities are based on the Uppsala package 
(Gustafsson 1973 1 which comprises continuous opacities only. 
Several Fe bound-free extinction edges though are located in 
the ultraviolet where a large number of atomic lines from met- 
als contribute to the UV "line haze" typically encountered in 
the spectra of cool stars. Accurate evaluation of the photo- 
ionisation rates in correspondence with these bound-free tran- 



2 loge Fc = log(j? F e/«H) + 12. 
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Table 2. Fe abundance corrections Alog e Fe = (log e Fe ) non -LTE _ (l°g ^Fe) lte f° r our sample stars. The results for calculations 
with and without line-blocking are shown, for two different values of the efficiency Su of Hi collisions. The reported values of 
Alog e Fe are averages over weak Fei lines in the visible and near UV with equivalent widths between 5 mA and 100 mA. The 
standard deviation is given as a measure of the line-to-line scatter in the non-LTE abundance corrections. No direct comparison 
with observations has here been undertaken. 





Fe 


non-LTE Abundance Corrections from Fe i Lines 






S 


h = 0.001 


S H = 1 




Model 


No blocking 


With blocking 


No blocking 


With blocking 


Sun (marcs) 


0.24 ± 0.02 


0.08 ± 0.02 


0.12 ±0.01 


0.03 ±0.01 


Sun (hm) 


0.22 ± 0.02 


0.07 ± 0.02 


0.13 ±0.02 


0.04 ± 0.01 


Sun (val-3c) 


0.25 ± 0.04 


0.12 ±0.03 


0.17 ±0.03 


0.07 ± 0.02 


Procyon 


0.23 ± 0.03 


0.09 ± 0.02 


0.13 ±0.02 


0.04 ± 0.02 


HD 140283 


0.58 ± 0.05 


0.56 ± 0.05 


0.32 ± 0.02 


0.30 ± 0.02 


G64-12 


0.69 ± 0.06 


0.68 ± 0.06 


0.48 ± 0.03 


0.47 ± 0.03 



sitions requires that we take into account the effects of line- 
blocking associated with the UV line haze. 

Some previous attempts have been made to include the 
effects of line-blocking in non-LTE calculations in multi. In 
their study of non-LTE formation of the solar spectrum of al- 
kali Bruls et al. ([1992 1 reproduced the contribution of atomic 
lines to the background opacities by multiplying the total H 
and continuous metal opacities with wavelength dependent 
factors. Their factors were determined empirically although 
only for the Sun by fitting the observed solar continuum with 
the emerging continuum intensity computed with the val-3c 
model solar atmosphere and assuming they are independent of 
depth. Andretta et al. ( 1997 1 introduced background atomic and 
molecular lines to the standard multi's opacities in the form of 
an opacity table provided with the particular model atmosphere 
they used. The opacity dependence upon wavelength was given 
in their table at each depth point of the model atmosphere with 
a resolution of 2 A. The opacity data they used were fixed 
however and model-specific. More recently Bayazitov (2003 1 
investigated the effects of line-blocking on non-LTE Fe i line 
formation by adding the contribution of about 1 8 000 spectral 
lines to multi's continuous opacities. His work however was re- 
stricted to stars with metallicities close to solar and relied on a 
much simpler model Fe i atom with 99 levels plus continuum. 

In our study, in order to include the effects of line-blocking 
in multi we proceeded by sampling metal line opacities for 
about 9 000 wavelength points in the spectral region between 
1 000 A and 20 000 A and adding them to the standard back- 
ground continuous opacities. In order not to count the contri- 
bution of Fe line opacities twice however, metal line opacities 
were only added to the background when computing intensi- 
ties for the photo-ionisation rates and not for Fe i and Fe n line 
profiles. Molecular line opacities were not considered in our 
line formation calculations which however is a reasonable ap- 
proach for our stars where the atomic contribution dominates 
overwhelmingly. We performed simple bench-mark test LTE 
radiative transfer calculations with background line opacities 
sampled between 9 000 and 70 000 wavelength points cover- 
ing the same wavelength region. The tests indicate that photo- 



ionisation rates computed with our choice of opacity sampling 
for 9 000 wavelengths are sufficiently accurate for our pur- 
poses. In addition, we performed non-LTE test calculations 
sampling the opacities at different wavelengths and found that 
the overall effect on the final Fe non-LTE abundance correc- 
tions is less then 0.01 dex. 

For the calculations presented in this paper the background 
metal line opacities were computed in detail by our modified 
version of multi, consistently with the particular model atmo- 
sphere given as input, as a function of temperature, gas- and 
electron-pressure at each depth point. For this purpose we used 
the opacity package from the os-marcs code (Gustafsson et al. 
2005 1 taking into account the specific chemical composition 
of the given model stellar atmospheres. Level populations and 
ionisation equilibria for the various metals - including Fe - 
were computed assuming LTE. This inconsistency for Fe is 
however of minor importance and may well be balanced by 
missing line-blocking for other species. 

The non-LTE abundance corrections were estimated by 
computing both the LTE and non-LTE cases for several dif- 
ferent abundances, from which the non-LTE Fe abundance that 
reproduced the LTE equivalent widths were interpolated. For 
each set of calculations with a given model atmosphere we se- 
lected the Fe i spectral lines in our model atom falling in the 
visible and in the near UV, (3 000 A < A < 10 000 A), and with 
equivalent widths W,\ between 5 mA and 100 mA. 

3. Results 

Figure|2]shows the non-LTE corrections to the iron abundance 
derived for the marcs model solar atmosphere for our selec- 
tion of Fe i lines; the behaviour for the other two solar models 
is similar as summarised in Table |2] Results are presented for 
the two different efficiencies of Drawinian collisions adopted 
in this study, Sh = 0.001 and Su = 1, and for both the case 
including line-blocking and the one excluding it. 

As a consequence of over-ionisation the overall line opacity 
of weak neutral iron lines decreases and the lines form deeper 
into the atmosphere in non-LTE than with the assumption of 
LTE, as already shown by Athay & Lites ( I1972> . The equiva- 
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Fig. 2. Non-LTE Fe abundance corrections for the Sun's (Top) and Procyon's (Bottom) marcs model atmospheres. The results 
are presented for two values for the efficiency S h of Hi collisions: Sh = 0.001 (Left) and Sh = 1 (Right). Empty circles refer to 
the calculations without background line opacities included, filled circles to the calculations with background line opacities. For 
the Sun the corrections computed including background line opacities a la Bruls et al. (1992) are also plotted (crosses). 



lent widths of these lines thus decrease and therefore a positive 
correction has to be applied to the LTE iron abundance derived 
from Fei lines. As expected, the non-LTE effects are more se- 
vere in models with a steeper temperature structure in the line 
forming layers as the val-3c model (see discussion in Asplund 
120051 . 

Let us first examine the abundance corrections in the case of 
low-efficiency Hi collisions. When background line opacities 
are not included in the calculations we estimate a correction 
to the iron abundance of about 0.24 dex for the weak (W^ 5 
50 mA) solar Fei lines in our sample with little dependency 
on the equivalent width W,\. Weak lines form generally deep 
in the stellar atmosphere where the source function S v is given 
by the local Planck function B Y (T) and the non-LTE equivalent 
width reflects mainly Fe i depopulation in the region of line 
formation. Stronger lines (W,\ <: 50 mA) on the contrary form 
further out where the source function at line centre typically 
falls below the Planck function compensating in part for the 
effects of Fei depletion on the line strength (Saxner 1984). 

With the inclusion of sampled background line opacities in 
the calculations the mean intensity field J v at a given depth de- 



creases and so do the radiative ionisation rates which are given 
by Rion - J ^n^Jvdv where <x v indicates the photo-ionisation 
cross-section. A lower ratio between radiative ionisations and 
recombinations and therefore a weaker over-ionisation are ex- 
pected. Consequently the non-LTE Fe abundance corrections 
also decrease and fall to values of about 0.08 dex (average over 
lines with W,\ between 5 and 100 mA). 

We investigated as well the effects of line-blocking when 
multiplication factors were applied to the continuum opacities 
to simulate the UV line haze according to the recipe of Bruls 
et al. ( 1992 1. The results indicate that over-ionisation is indeed 
reduced in this case but by a smaller amount: the predicted Fe i 
non-LTE abundance corrections remain about 0.05 dex higher 
than in our calculations with line-blocking. The reason for the 
difference is that the technique by Bruls et al. ( 1992i results in 
overall lower mean background opacities than in the case of the 
full opacity sampling method. 

When the rates of collisions with neutral hydrogen are in- 
creased (Sh = 1) a lower degree of over-ionisation is achieved 
as expected due to increased thermalisation. The equivalent 
widths of non-LTE weak neutral iron lines therefore increase 
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Fig. 3. Non-LTE Fe abundance corrections for the metal-poor stars HD 140283's (Top) and G64-12's (Bottom) marcs model 
atmospheres. The results are presented for two values for the efficiency Sh of Hi collisions: Sh = 0.001 (Left) and Sh = 1 
(Right). Empty circles refer to the calculations without background line opacities included, filled circles to the calculations with 
background line opacities. 



and the overall non-LTE Fe abundance correction is reduced as 
illustrated in Fig.|5]and in Table|2 Because of the overall lower 
non-LTE abundance corrections, the impact of background line 
opacities is less pronounced with a difference of about 0. 1 dex 
between the calculations including and excluding line-blocking 
at solar metallicity. 

Concerning Procyon (Fig. O, departures of Fe abundance 
from LTE calculations are similar to the solar case with 
and without line-blocking (Table |3. Because of its higher 
r e jf and lower logg Procyon shows slightly larger non-LTE 
iron abundance corrections than the Sun for weak Fei lines. 
Multiplication factors from Bruls et al. ( 1992 1 were calibrated 
specifically for the solar atmosphere and therefore were not ap- 
plied to simulate line-blocking in Procyon's atmosphere. 

In contrast to the solar metallicity case, the inclusion of 
the line haze at low metallicities produces minor changes in 
terms of mean intensity J v due to the overall low contribution 
of lines to the total opacities. The exclusion of background line 
opacities still produces systematically larger mean intensities 
at all depth and therefore generally increases photo-ionisation 
rates and non-LTE Fe abundance corrections but only very 



slightly so. The difference between the results of calculations 
with and without line-blocking is typically only 0.02 dex for 
both HD 140283 and G64-12. 

Figure @] illustrates the effects of line-blocking on the de- 
parture from LTE of total Fei number density as a function 
of depth for the marcs model atmospheres of the Sun and the 
metal-poor HD 140283; the departure coefficients for individ- 
ual Fe i levels are all similar and considerable in the latter case. 
At solar metallicity the inclusion of line-blocking significantly 
reduces the departure of the Fe i population from LTE in the re- 
gion where most weak lines form. On the contrary at low metal- 
licity calculations with and without line-blocking produce very 
similar and significant non-LTE signatures for a large portion 
of the line forming region in the atmosphere. Due to the overall 
low contribution of the line haze to the opacities at low metal- 
licities, the mean intensities are typically larger than at solar 
metallicity. Consequently photo-ionisation rates are higher and 
departures from LTE more pronounced in metal-poor stars. In 
fact the departure of Fe i population from LTE in metal-poor 
stars is controlled primarily by the efficiency of the H i colli- 
sions. Changing the scaling factor S h in Drawin's formula from 
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Table 3. Fe non-LTE abundance corrections for our sample stars estimated from weak Fe n lines in the visible and near UV. 
Averages are computed over weak Fe n lines in the visible and near UV with equivalent widths between 5 mA and 100 mA. The 
standard deviation is given as a measure of the line-to-line scatter in the non-LTE abundance corrections. 



Fe non-LTE Abundance Corrections from Fe n Lines 
5 H = 0.001 5 H = 1 



Model 


No blocking 


With blocking 


No blocking 


With blocking 


Sun (marcs) 


0.02 ± 0.01 


0.01 ± 0.01 


0.00 ± 0.00 


0.00 ± 0.00 


Sun (hm) 


0.02 ± 0.01 


0.01 ±0.01 


0.00 ± 0.00 


0.00 ± 0.00 


Sun (val-3c) 


0.00 ± 0.00 


0.00 ± 0.00 


0.00 ± 0.00 


0.00 ± 0.00 


Procyon 


0.02 ± 0.01 


0.02 ± 0.01 


0.01 ±0.00 


0.01 ± 0.00 


HD 140283 


0.24 ±0.13 


0.24 ±0.13 


0.01 ±0.01 


0.01 ±0.01 


G64-12 


0.17 ±0.08 


0.17 ±0.08 


0.05 ± 0.06 


0.05 ± 0.06 




Log(T 5D00 ) Log(r 5000 ) 

Fig. 4. Departure of total Fei number density rip eI from LTE versus standard optical depth at A = 5000 A for the Sun (Left) 
and for HD 140283 (Right). non-LTE departures are shown for the case of fully efficient Drawinian collisions and for different 
treatments of line-blocking. 



0.001 to 1.0 decreases the non-LTE abundance corrections by 
0.2 - 0.3 dex for the metal-poor stars studied here. 

It is clear that taking account of line-blocking is impor- 
tant for the Fe i non-LTE line formation at solar metallicity but 
significantly less so at low metallicity. Instead, the dominant 
source of error in metal-poor stars is the treatment of inelastic 
H i collisions, with the non-LTE abundance corrections differ- 
ing by 0.2 - 0.3 dex when using Sh = 1 and Sh — 0.001 (Table 
|2j. While our results for the metal-poor case with Sh — 0.001 
are very similar to the findings of Korn et al. ( 2003 1 with no hy- 
drogen collisions, we note that our Fe non-LTE abundance cor- 
rections with Sh — 1 are still significantly larger than in their 
calculations. The main reason for this is that in addition to us- 
ing even more efficient Hi collisions (Sh = 3), they also more 
importantly enforce ad-hoc thermalisation of the uppermost 
levels of Fe i with the ground level of Fe n. Since the latter is the 
dominant Fe state and is well described by LTE, all Fe i levels 
are steered closer to the LTE solution by bound-bound colli- 
sions within the Fe i system. In fact test calculations performed 
with their model atom for HD 140283 and with Sh = 1 but no 
enforced thermalisation of highly excited Fe i levels produce 
similarly large non-LTE effects, i.e. Alog e Fe « 0.45 dex com- 



pared with our value of 0.3 dex (Korn 2005, private commu- 
nication). Our non-LTE abundance corrections are also much 
larger than the corresponding computations of Gratton et al. 
( 1999 1, which is not surprising given their very efficient H i col- 
lisions (Sh ~ 30) and their use of considerably smaller photo- 
ionisation cross-sections than those of Bautista et al. ( 1997 1. 

Overall we find slightly larger non-LTE corrections than 
Thevenin & Idiart ( 1999 1 even when we consider fully efficient 
H i collisions (Sh = 1) as it was actually used in their work. On 
the other hand the approach of Thevenin & Idiart ( 1999 1 is not 
identical to ours. Their Fei and Fen systems were not merged 
and their analysis relied on a different suite of model atmo- 
spheres. As mentioned before our treatment of line-blocking 
differs from theirs as they used multiplication factors to ac- 
count for the background line haze in a similar way as Bruls 
et al. ( 1992 1 did for the solar case. 

Finally concerning Fe n lines, Table [5] shows the Fe non- 
LTE abundance corrections estimated for weak singly ionised 
iron lines in the visible and near UV. Contrary to what is of- 
ten assumed Fen lines are not exempt from non-LTE effects. 
Departures of Fe n lines from LTE can be substantial especially 
at low metallicity. In particular, in the case of low H i colli- 
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sions Fe non-LTE abundance corrections of about 0.2 dex are 
found for metal-poor stars. Similar corrections for HD 140283 
are also reported by Shchukina et al. (2005 1. For fully efficient 
H i collisions corrections are however much lower although not 
completely vanishing at low metallicity. We note that the net 
Fe non-LTE effect (i.e. Aloge FeI - Aloge FeII ) is quite similar 
for S h — 0.00 1 and S h = 1 . It may therefore be difficult to em- 
pirically discriminate between different choices of S h on the 
basis of observations alone. 

4. Conclusions 

We have investigated the effects of line-blocking on non-LTE 
Fe abundances derived from Fe i lines for a selection of four 
representative late-type stars. We found the effects of line- 
blocking to be significant at solar metallicities. In particular 
calculations including background line opacities result in non- 
LTE Fe abundance corrections 0.1 - 0.15 dex higher than in 
calculations excluding them, depending on the efficiency of Hi 
collisions a la Drawin (S^ = 1 and Su = 0.001 respectively). 
On the other hand line-blocking has a much smaller impact at 
low metallicities. Absolute non-LTE Fe abundance corrections 
are in this case significant and sensitive mainly to the strength 
of H i collisions and vary only by 0.01 - 0.02 dex between cal- 
culations with and without line-blocking. At the present time 
the main uncertainty in non-LTE line formation calculations in 
metal-poor stars is still the treatment of the poorly known in- 
elastic H i collisions and whether or not thermalisation of the 
highly excited Fe i levels should be applied. 
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